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Abstract 

The atmospheric neutrino flux calculated with FLUKA was originally limited 
to 100-^200 GeV for statistical reasons. In order to make it available for the analysis 
of high energy events, like upward through-going muons detected by v telescopes, 
we have extended the calculation so to provide a reliable neutrino yield per primary 
nucleon up to about 10 6 GeV/nucleon, as far as the interaction model is concerned. 
We point out that the primary flux model above 100 GeV/nucleon still contributes 
with an important systematic error to the v flux. 

1 Introduction 

The main motivation of FLUKA[1] calculation of atmospheric neutrino flux[2] is the 
attempt to minimize as much as possible the theoretical uncertainties connected to the 
shower and hadronic interaction models. The high accuracy of the algorithms adopted 
in FLUKA and the wide set of experimental data used to benchmark its hadronic and 
electro-magnetic interaction models have already allowed the use of this simulation mod- 
els in many applications where the maximum available precision is requested. FLUKA 
has been successfully used also in cosmic ray physics. For this reason it is also going 
to be adopted as an option in the CORSIKA code, at least for interactions below 80 
GeV[3]. Our results concerning the atmospheric neutrino flux calculations have been pre- 
sented in [2] and references therein . There, the other important factors contributing to 
the systematic uncertainties which are not connected to the FLUKA code, are the atmo- 
sphere description, the accuracy of the geomagnetic model, and, as the most important 
contribute, the knowledge of the primary flux. 

The practical problems arising from our approach is that the demand of computer 
power for very high energy showers may be very large. This is the reason why so far we 



have presented calculation results only up to about 100 GeV (useful for the analysis of 
contained events) due to insufficient statistics at higher energy. The low energy sector 
was also considered less affected by the uncertainty in the primary flux thanks to the high 
accuracy data from AMS[4] and BESS [5]. This statistical limitation in energy has however 
prevented the use of FLUKA results for a complete analysis of atmospheric neutrinos, 
including through-going muons produced by neutrinos with mean energies 50-100 GeV. 
It has been already shown how this requires the knowledge of neutrino flux up to at least 
10 4 GeV. In the high energy sector the existing available calculations are those coming 
mostly from faster approaches, like those of Bartol[6] and HKKM[7]. Also the new recent 
HKKM calculation [8] (HKKM2001) stops in practice at about 1 TeV and makes use of 
extrapolations and/or simplifications. 

In order to overcome this limitations, we have performed an extensive simulation run 
so to extend the FLUKA calculation of atmospheric neutrino fluxes up to 10 4 GeV. This 
has been done in such a way to consider different possible primary spectra. Results are 
presented in this paper. 

2 The Calculation Model 

We have used the same calculation set-up reported in [2], where also some details about 
the physical model are reported. Here the different strategy is to generate neutrinos 
avoiding the sampling from a unique primary all-nucleon spectrum and instead prepare 
parallel simulation runs for many different small energy intervals, summing together the 
results with a proper weighting scheme so to reproduce any desired input spectrum. In 
this way we could obtain a statistical uncertainty almost independent on energy. We 
started from few tens of GeV of primary energy and in order to achieve reliable results at 
a maximum neutrino energy of 10 4 GeV, we have pushed the maximum nucleon energy up 
to 10 6 GeV. The relatively high energy allowed some simplification, avoiding the question 
of solar modulation and geomagnetic cutoff. Furthermore, the angular distribution for 
non oscillated flux has been recorded in 10 bins in a single hemisphere. 

3 Results and Discussion 

We have chosen two different primary all-nucleon spectra to present and compare the 
calculation results. The first one is the 2001 Bartol fit of ref. [9] taking into account the 
most recent direct measurements, and that was used as reference in our last work[2]. The 
second one is the previous Bartol fit to all-nucleon flux of ref. [6] , which was already used 
by us for the first (before 2001) FLUKA calculations. The updated flux tables, obtained 
using the 2001 Bartol fit, are available from our web site [10]. Numerical values of the 
and flux above 20 GeV are reported in Tables 1 and 2. In Tables3 and 4 we instead 
report the neutrino flux obtained using the old Bartol fit to primary flux of ref. [6]. 

In Fig.l we show the angle integrated flux for muon v and V for the two primary 
spectra. It is evident that in the high energy sector the choice between the two different 
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options for the primary all-nucleon spectrum has a large impact. 

In order to give a quantitative figure of this statement in Fig. 2 we also show, as a 
function of neutrino energy, the ratio of the flux calculated using the old Bartol fit with 
respect to the one obtained with the 2001 fit (solid line) for and v^. In the same plots 
we also show the ratio of the Bartol and HKKM2001 flux with respect to FLUKA with 
the 2001 primary fit. 

The difference introduced by using the two primary spectra evidently results from 
a different weight given to the high energy primary measurements. The highest energy 
point in the recent high precision experimental measurements from AMS and BESS is 
at 100 GeV. Higher energy points in the TeV region (essentially from JACEE[11] and 
RUNJOB[12]) have much larger error bars. There are indeed suspicions that the multi- 
TeV proton component could be underestimated in the 2001 fit. Assuming to consider 
the 100 GeV point as a pivot from which a simple power law spectrum (-E"" 7 ) is used to 
extrapolate to the high energy region, we obtain the results shown in Fig.3, by varying 
the spectral index 7 from 2.6 to 3.0. We also show the flux ratios. 

4 Conclusions 

The FLUKA atmospheric neutrino spectrum calculation has been extended up to 10 4 
GeV, so that it can be used for the analysis of upward going muon experiments. Above 1 
TeV, the uncertainty in the primary flux seems to be the most important contribution: it 
may exceed 20%. In principle this has no impact on the essential feature of experimental 
analysis, namely the deformation of the angular distribution of upward going muons as 
induced by oscillations. However, also the absolute normalization has some importance, 
at least to achieve a better internal consistency. As a first guess, considering the event rate 
of upward-going muons observed by MACRO [13] and Super-Kamiokande[14], the results 
that we obtain using the 2001 Bartol fit seems to exhibit a too low normalization (of the 
order of 15-i-20 %, if we include oscillations), although this has to be confirmed by detailed 
analysis performed by the experimental collaborations. Within the statistical uncertainty 
of the 2001 primary fit, our results are very close to those of the HKKM2001 calculation. 
A less steep primary spectrum at high energy, as that of ref. [6], would probably provide 
a better absolute normalization. 
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Figure 1: Angle integrated (and without oscillations) atmospheric and FLUKA 
fluxes for the 2 choices of primary spectrum considered here, weighted by E 2,5 . 
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Figure 2: Left: flux ratio vs. energy with respect to the FLUKA + 2001 primary fit 
of: 1) the FLUKA+Old Bartol primary fit (solid line), Bartol flux (dotted green line) and 
HKKM2001 flux (dashed red line). Right: the same for v^. 
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Figure 3: Left: Angle averaged FLUKA fluxes weighted by E 3 , including also the 
case of simple power law primary spectra with different spectral indexes, ranging from 2.6 
(upper curve) to 3.0 (lower curve). The HKKM2001 flux is also shown. Right: ratios of 
the different cases with respect to the present FLUKA result. 
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Table 1: Muon neutrino flux in units of (cm 2 s sr GeV) -1 above 20 GeV as a function of 
energy for different values of the cosine of zenith angle, for the case of the 2001 Bartol fit 
of primaries [9]. 
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Table 2: Muon anti- neutrino flux in units of (cm 2 s sr GeV) -1 above 20 GeV as a function 
of energy for different values of the cosine of zenith angle, for the case of the 2001 Bartol 
fit of primaries [9]. 
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Table 3: Muon neutrino flux in units of (cm 2 s sr GeV) -1 above 20 GeV as a function of 
energy for different values of the cosine of zenith angle, for the case of the old Bartol fit 
of primaries [6]. 
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Table 4: Muon anti- neutrino flux in units of (cm 2 s sr GeV) -1 above 20 GeV as a function 
of energy for different values of the cosine of zenith angle, for the case of the old Bartol 
fit of primaries [6]. 
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